Temporally heterogeneous dynamics in granular flows 
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Granular simulations are used to probe the particle scale dynamics at short, intermediate, and 
long time scales for gravity driven, dense granular flows down an inclined plane. On approach to the 
angle of repose, where motion ceases, the dynamics become intermittent over intermediate times, 
with strong temporal correlations between particle motions - temporally heterogeneous dynamics. 
This intermittency is characterised through large scale structural events whereby the contact network 
periodically spans the system. A characteristic time scale associated with these processes increases 
as the stopped state is approached. These features are discussed in the context of the dynamics of 
supercooled liquids near the glass transition. 

PACS numbers: 81.05.Rm,45.70.Mg,64.70.Pf,61.20Lc 



Granular materials persist at the forefront of contem- 
porary research due to the extreme richness in their dy- 
namics, either under shear 0], flowing out of a hopper 0, 
or driven by gravity Q. Because thermodynamic tem- 
perature plays little role in determining these features, 
granular materials are widely recognized as a macroscopic 
analogue of athermal systems far from equilibrium |^. 
The notion of jamming provides a unifying frame- 
work to describe the behaviour of a wide range of sys- 
tems, from the molecular level (supercooled liquids and 
glasses), microscale (colloidal glasses), through to the 
macroscale (granular matter), near their jamming transi- 
tion. The motivation for this work comes from the view 
that granular flows can be used as a macroscopically ac- 
cessible analogue through which we can gain a better 
understanding of the critical slowing down of the dynam- 
ics in systems near their point of jamming. One of the 
features associated with the approach of the glass tran- 
sition is the concept of spatially heterogeneous dynamics 
. This work presents features of a driven granular ma- 
terial flowing down an inclined plane - chute flows - that 
exhibit large scale, cooperative events on approach to the 
angle of repose where flow stops - its jamming point. 

Gravity driven flows continue to be studied extensively 
due to their ubiquity throughout nature - avalanches, de- 
bris flows, and sand dunes - and also because of their 
importance in materials handling throughout industry. 
More recently, there has been a focus of attention to- 
wards understanding the properties of chute flows be- 
cause of the ability of obtaining well-defined and repro- 
ducible steady state flow regimes. Despite the fact that 
recent experiments have been able to determine, in detail, 
macroscopic information of granular flows 0, 1^ , probing 
the internal dynamics remains problematic due to the 
opaqueness of the particles. Still, such flows offer a rela- 
tively clean system in which to develop theoretical models 
for constitutive relations as well as describing the flow at 
the level of the grain size j^. This is where the role of 
simulation has proved extremely useful. 



Granular simulations of chute flows have captured fea- 
tures of the rheology with remarkable accuracy [lol . 
Three principal flow regimes are observed for chutes with 
rough bases (see Fig.^l. If 6 is the inclination angle of the 
chute: i) no flow occurs below the angle of repose 6r, ii) 
steady state flow exists in the region 9r < < 9u, and iii) 
unstable flow occurs for 9 > 9u- For all 9 < 9u, the den- 
sity is constant throughout the bulk, away from the free 
top surface and lower boundary, with smoothly varying 
depth profiles of the flow velocity and the principle and 
shear stresses JJi] . The location of 9r also depends on the 
height of the flowing pile h 0, 0| ■ Pouliquen introduced 
the quantity hstop{0) that encodes all the (undetermined) 
information about the roughness and effective friction of 
the base to relate the dependency of 9r on h. In the 
language of glasses and colloids, hstopi^) represents the 
glass transition or yield stress line, respectively. 




FIG. 1: Phase behaviour of granular chute flow (schematic). 
9r is the angle of repose at which flow ceases. 9m is the maxi- 
mum angle of stability at which flow is initiated from rest. 9u 
is the maximum angle for which stable flow is observed. The 
angle 9i at which intermittent dynamics emerges, the shaded 
region, coincides with 9i k, 6m- For the simulations reported 
here, 9r = 19.30° ± 0.01°, 9m = 20.3° ± 0.2°, and 9^ ^ 26°. 



This work characterises the large-scale fluctuations and 
emerging spatial and temporally heterogeneous dynam- 
ics for a flowing state in the vicinity of 9^- This is in 
contrast to discrete avalanches in heap flows 14] . 9i is 
the angle, for 9r < 9 < 9i, such that, over time scales in- 
termediate between ballistic propagation, at short times, 
and diffusive motion, at long times, the particle dynam- 
ics become strongly intermittent. This intermittency is a 
consequence of large spatial and temporal correlations in 
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particle motions, through the relaxation and reformation 
of mechanically stable clusters that periodically span the 
system. A characteristic time scale r^, associated with 
these correlations is found to increase as [9 — 9^) — > 0. 

The simulations were performed using the discrete ele- 
ment technique with interactions appropriate for granu- 
lar materials Monodisperse (soft) spheres of diame- 
ter d and mass m flow down a roughened, inclined plane 
due to gravity g. Initially, 9 is set to 9u to remove any 
preparation history, then incrementally reduced in steps 
of 0.01° < A6' < 0.5°. Particles interact only on con- 
tact (cohesionless), along directions normal and tangen- 
tial (via static friction) to the vector joining their centres. 
Contacts are defined as two overlapping neighbours. All 
runs were performed with a particle friction coefficient 
/i = 0.5 and coefficient of restitution e = 0.88. Most of 
the results presented are for N — 8000 particles flowing 
down a chute that is = 2Qd long, Ly ~ IQd wide (the 
xj/-plane employs periodic boundary conditions), and an 
average flow height h « AOd (with a free top surface). See 
Fig. Massively, large-scale, parallel simulations were 
used for N = 160000, with {L^,Ly,h) = (100, 40, 40)d, 
to demonstrate that the observed features are not a sys- 
tem size artifact. Length scales are non-dimensionalised 
by the particle diameter d and time in units of t = yfdfg. 

The intermittent regime was first identified through 
measurements of the kinetic energy as shown in Fig. ^ 
In the continuous-flow, steady state regime the energy 
is approximately constant with small fluctuations about 
the mean value. As the inclination angle is incrementally 
decreased towards 9r^ changes in the time profile of the 
energy first become apparent for 9i k. 20.2°. 
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FIG. 2: Average kinetic energy per particle over an intermedi- 
ate time window for A*' = 8000 and different inclination angles 
as indicated. The symbols are for iV = 160000 at S = 19.5°. 

Because of the convective nature of the flow, the mean- 
squared displacement (MSD), shown in Fig. 13 is mea- 



sured normal to the inclined plane in the z direction: 
Az^ =< (^(t) — 2;(0))^ >. During continuous flow, 
there is a clear crossover between the ballistic regime, for 
t/r < 10~^, and diffusive regime, tjr > 10. Closer to the 
jamming point, the MSD oscillates over times intermedi- 
ate between them. The appearance of novel dynamics in 




FIG. 3: The mean squared displacement in the direction nor- 
mal to the plane A^;^, with decreasing angle from top to bot- 
tom. All curves have been shifted for clarity. 

the vicinity of 9t is a consequence of strong correlations 
emerging in time-time quantities such as the velocity au- 
tocorrelation function, C{€) =< Vz{t)vz{0) >, where Vz 
is the velocity component in the direction normal to the 
plane. See Fig. 2Ib). These fluctuating quantities are a 
signature of temporally heterogeneous dynamics. To in- 
vestigate this further, the average coordination number 
z, Fig. 01a), is shown over the intermediate time win- 
dow, demonstrating that the structural processes asso- 
ciated with these correlations involve large-scale cooper- 
ative events. During the flowing phase of the intermit- 
tent regime, z remains lower than the mechanically stable 
limit for frictional spheres Zc = 4 [isj . In the static phase, 
the systems attains a value z m Zc, thus almost coming 
to rest. (There is always a residual creeping motion as 
flow does not completely stop until 9r-) 

To facilitate such processes, the contact force network 
undergoes temporal transitions between static and flow- 
ing phases. This evolution is captured in the (distinct) 
particle-particle contact force-force time correlation func- 
tion J^(i), as shown in Fig^Ic). The physical picture of 
the intermittent regime is thus: over some characteris- 
tic time scale the system oscillates between almost total 
mechanical stability and fluid-like flow. This involves 
time-dependent, system-spanning, structural relaxation 
events through the break-up and reformation of particle 
contacts. This process is best depicted in the simulation 
snapshots of Fig. \S\ 

The intermittent regime can therefore be considered as 
temporally biphasic. During one phase, z « z^, so the 
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FIG. 4: Intermediate time window associated with the oscilla- 
tions in the MSD, for 19.35° (dashed line), 19.5° (thick), and 
24° (dotted): (a) coordination number z, (b) velocity auto- 
correlation function C(€), and (c) force- force time correlation 
function T(t). 



depending on the time window of observation. 
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FIG. 6: Probability distribution function P(/) of the normal 
contact forces (normalised so / = 1) for Q = 20°; total (solid 
line), the high-z phase (dashed), and the low-z phase (dot- 
ted). In the high-z phase, the contact forces are binned into 
a histogram if 2: > 3.5 for that configuration, and similarly 
\i z < 2.0 in the low-2 curve. The full P(/) for Q = 22° 
(symbols) is shown for comparison. 




FIG. 5: The contact force network at S = 20° for three suc- 
cessive times, (a) Snapshot corresponds to a peak in the co- 
ordination number, (b) a trough, and (c) the next peak. The 
shade and thickness of the lines represent the magnitudes of 
the normal forces between two particles in contact. The in- 
clination has been removed in the figures and the frame is a 
guide to the eye. The simulation axes are indicated. 



system is close to mechanical stability, while the flowing 
phase has a considerably lower coordination. A signature 
of this is seen in the probability distribution function of 
the normal contact forces P{f)- The partial distribu- 
tions shown in Fig. are for the high-z phase (defined 
as configurations with z > 3.5) which resembles that of 
a static packing - exponential decay at high forces - and 
separately the low-z phase (z < 2), with a high-/ tail 
that decays more slowly. This suggests that it may be 
possible to distinguish between different types of flows 



The intermittent regime is peculiar due to the complex 
interplay between several different characteristic length 
and time scales. Identifying the structural entities of 
characteristic size £, remains a matter of debate 0, ■ 
However, the general picture is as follows: one expects 
£ to increase with decreasing 9, or correspondingly, the 
shear rate 7. Then, 9^. is the angle where i h, resulting 
in the transition from a dynamic state to a static one. 
The development of an intermittent regime suggests an- 
other competing time scale, whereby £ ^ h can occur, 
but other mechanisms drive the system away from sta- 
bility before the system permanently jams. This seems to 
coincide with the regime where the scaling of the flow ve- 
locity differs from that of the continuous flow regime ■ 
Presently, it is not clear how to approach this problem 
without a more sophisticated particle scale analysis. 

Progress can still be made by identifying the charac- 
teristic time scales of the intermittent regime. One is the 
periodicity in the kinetic energy, correlation functions, 
and z. This time scale is rather insensitive to 9 in the in- 
termittent phase. However, C{t) has an envelope of decay 
that corresponds to the extent of the intermittency. A 
time scale Tc, that characterises the initial decay in C{t) 
is shown in Fig[7| Interestingly, a power-law divergence: 
Tc ^ (0 — 0r) with 7 « 0.65, captures the behaviour 
quite well. As a comparison, a modifled Vogel-Fulcher 
form is also shown in Fig.[7| which is satisfactory close to 
the jamming point, but deviates further away. The avail- 
able range for analysis is somewhat restricted through 
this particular deflnition of Tc, so it remains unclear how 
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this time scale might correspond to the more famiHar def- 
initions of relaxation times in supercooled fluids JJ^. A 
more appropriate definition of is being investigated. 




FIG. 7: Characteristic time scale Tc, with error bars, asso- 
ciated with the intermittent structural events at intermedi- 
ate times as a function of distance from the jamming point 
{9 — 6r). The solid line is a power law fit: t ^ (6 — Or) 
with 7 ~ 0.65; whereas the dashed line is of the modified 
Vogel-Fulcher form: Tc ~ exp{A{e - 6'r)"''}, with A « 0.66 
and p ~ 0.39. 

In conclusion, gravity-driven, dense, granular flows 
down an inclined plane exhibit intermittent dynamics 
over intermediate time scales in the vicinity of the an- 
gle of repose. The intermittent regime signals the onset 
of temporal heterogeneous dynamics, whereby the short 
time motion is ballistic, at long times diffusive, but at 
intermediate times a combination of creeping flow and 
trap-like dynamics exists. The system becomes trapped 
in temporally metastable states jlq through large scale 
cooperative events. Correlations between these events 
decay, but do so more slowly as the jamming point is ap- 
proached. The behaviour of the characteristic time scale 
Tc, shares similarities with the properties of supercooled 
liquids in the vicinity of the glass transition. Such dy- 
namical similarities were recently probed experimentally 
in vibrated bead packs ji^. Therefore, in a typically 
heuristic fashion, an analogy with glassy dynamics the- 
ory can be made. Or marks the point where all dynamics 
cease and in this sense plays the role of the glas s transi- 
tion temperature Tg. Indeed, earlier work showed 
that the flowing-static transition resembles the way a 
model liquid approaches the glassy phase 0. However, 
the way the dynamics change at 9i , is similar is spirit to 
the way the dynamics in a supercooled fluid changes at 
the mode coupling temperature Tc 0| . With this view in 
mind, one can associate 9i with Tc- These results there- 
fore hint towards a unifying picture of dynamical het- 
erogeneities in dense, amorphous systems. Because the 



time scales associated with the onset of dynamical het- 
erogeneity and cooperativity are experimentally accessi- 
ble, granular flows offer a convenient system to address 
questions often associated with glasses at the molecular 
and colloidal level. This emerging picture thus begs the 
question, can existing theories that are currently being 
applied to traditional glasses and colloids be extended to 
granular systems? 
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